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We investigate the implications of models that achieve a Standard Model-like Higgs 
boson of mass near 125 GeV by introducing additional TeV-scale supermultiplets in 
the vector-like 10+10 representation of SU(5), within the context of gauge-mediated 
supersymmetry breaking. We study the resulting mass spectrum of superpartners, 
Q_i. comparing and contrasting to the usual gauge-mediated and CMSSM scenarios, and 

discuss implications for LHC supersymmetry searches. This approach implies that 
exotic vector- like fermions t\ 2 , b' , and r' should be within the reach of the LHC. We 
discuss the masses, the couplings to electroweak bosons, and the decay branching 
ratios of the exotic fermions, with and without various unification assumptions for 
the mass and mixing parameters. We comment on LHC prospects for discovery 

■ of the exotic fermion states, both for decays that are prompt and non-prompt on 

■ detector-crossing time scales. 
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I. INTRODUCTION 

Recently the ATLAS and CMS experiments at the LHC haveput forward data analysis results 
that suggest the Higgs boson mass could be close to 125 GeV [H, E|- The statistical significance 
is not at the 'discovery level', nor is it enough to determine if the putative Higgs boson signal is 
really that of the Standard Model (SM) Higgs boson, or some close cousin that may have somewhat 
different couplings and rates. Nevertheless, we wish to investigate the supposition that the Higgs 
boson exists at this mass and is SM-like in its couplings. 

Stipulating the above, supersymmetry is an ideal theoretical framework to accommodate the 
results. The many favorable features of supersymmetry are well-known [jj, but the one most appli- 
cable here is its generic prediction for a SM-like Higgs boson with mass less than about 130 GeV. 
Within some frameworks of supersymmetry, such as 'natural' versions of minimal supergravity 
(mSUGRA) or minimal gauge mediated supersymmetry breaking (GMSB), a Higgs mass value 
of '-mGeVt seems perhaps uncomfortably high. Within other frameworks, such as 'unnatural' 
PeV-scale supersymmetry [4!] or split supersymmetry [5J, such a mass value seems perhaps uncom- 
fortably low. Nevertheless, almost any approach to supersymmetry allows one to easily absorb this 
Higgs mass into the list of defining data and then present the resulting allowed parameter space. 

In this article we wish to see how well one can explain a ~ 125 GeV Higgs boson mass using 
'natural' supersymmetry. There are many good discussions of this already present in the liter- 
ature [{], 0] , but the approach we take here is to use extra vector-like matter supermultiplets to 
raise the Higgs mass [8(- 2lJ. As shown in detail in 12], the Yukawa coupling of the vector- like 



quarks to the Higgs has a fixed point at a value large enough to substantially increase the lightest 
Higgs mass while giving a fit to precision electroweak oblique observables that is as good as, or 
slightly better than, the SM. This can be done in various different scenarios for the soft terms, but 
here we choose to investigate within the context of GMSB; earlier studies of this can be found in 
The details of the specific model we study will be discussed in the next section. We 



like this approach because the superpartner masses are not required to become extremely heavy to 
raise the light Higgs mass through large logarithms in the radiative corrections, nor does one need 
to invoke very large Higgs-stop-antistop supersymmetry-breaking couplings. Instead, the extra 
vector-like states, interacting with the Higgs boson, make extra contributions to the Higgs boson 



mass in a natural way. This approach has been reemphasized recently also by [13]-[2l|] within the 
context of the ~ 125 GeV Higgs boson signal, and our study confirms some previous results and 
extends the understanding by investigating correlations within a unified theory and detailing the 
phenomenological implications that can be useful for the LHC experiments to confirm or reject 
this hypothesized explanation for the Higgs boson mass value. 



II. MINIMAL GMSB MODEL WITH EXTRA VECTOR-LIKE PARTICLES 



A. Theory definition, parameters and spectrum 

The theory under consideration here is a minimal GMSB theory with one SU (5) 5 + 5 messenger 
multiplet pair, along with a 10 + 10 multiplet pair at the TeV scale. We choose this model because 



In this article, ~125 GeV always means any value calculated theoretically to be between 122 GeV < Mh < 128 GeV, 
consistent with the LHC results taking into account experimental uncertainty (notably the lack of a definitive 
signal) as well as theoretical errors in calculating the Higgs mass from the supersymmetric input parameters. 
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it is minimal, illustrates the key phenomenological features of this broad class of theories, and 
maintains perturbative gauge coupling unification at the high scale. The same model has also been 
considered in 19M21I]. The unification scale (defined as the scale where g\ and g% meet) turns 



out to be larger than the corresponding scale in the MSSM by a factor of 2-4, depending on the 



sparticle thresholds and the GMSB messenger scale. As in ref. [I2j, we use 3-loop beta functions 
for the gauge couplings and gaugino masses, and 2-loop beta functions for all other parameters. 
These renormalization group equations are not given explicitly here, because they can be obtained 
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in a straightforward and automated way from the general results given in refs. 

To set the notation, the MSSM fields are defined below along with their SU(3)c x SU{2)l x 
U(1)y quantum numbers: 



qi = (3,2,1/6), ^ = (3,1,-2/3), di = (3,1,1/3), 

£i = (1,2,-1/2), ei = (1,1,1), fl„ = (1,2, 1/2), H d = (1, 2, -1/2). (2.1) 

with i = 1,2,3 denoting the three families. The MSSM superpotential, in the approximation that 
only third-family Yukawa couplings are included, is 

Wmssm = ^H u H d + y t H u q 3 u 3 - ybH d q 3 d 3 - y T H d £ 3 e 3 . (2.2) 

The 10 and 10 SU(5) multiplets are comprised of Q, U, E and Q, U, E supermultiplets, 
respectively, with 

Q = (3,2,1/6), U = (3,1,2/3), E = (1,1,-1), (2.3) 
Q = (3,2,-1/6), U = (3, 1,-2/3), E= (1,1,1). (2.4) 

These extra fields interact with the MSSM Higgs bosons at the renormalizable level. The relevant 
superpotential is 

W QUE = M Q QQ + MuUU + M E EE + kH u QU - k'H d QU. (2.5) 

The extra superfields of the 10 + 10 give rise to additional exotic particles beyond the MSSM: 
charge +2/3 quarks t' 12 (plus scalar superpartners ^1234), a charge —1/3 quark b' (plus scalar 
superpartners b\ 9), and a charged lepton r' (plus scalar superpartners f{ 2 )- 



As noted in 12j, the Yukawa interaction k is subject to an infrared-stable quasi-fixed point 
251 ] slightly above k = 1.0 at the TeV scale. This value is both natural (since a large range of 
high-scale input values closely approaches it), and is easily large enough to mediate a correction 
to the lightest Higgs boson mass that can accommodate ~125 GeV or larger, depending of 
course on the other parameters of the theory. In this paper, we will always assume that k is near 
its (strongly attractive) quasi-fixed point, by arbitrarily taking k = 1 near the apparent scale of 
gauge coupling unification and evolving it down. Taking larger values at the high scale would only 
increase the TeV-scale value of k by about 2% at most, although it should be kept in mind that 
the contribution to the Higgs squared mass correction scales like fc 4 . For simplicity, we will take 
k! to be small, since it does not help to raise the h° mass, although a small non-zero value would 
not affect the results below very much. The superpartner spectrum of this theory is determined by 
the normal procedures for minimal GMSB. The input parameters needed are tan/3, sign(/x), the 
mass scale for the 5 + 5 messenger masses M mess and the supersymmetry breaking transmission 
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scale A which is equal to (Fs)/(S) where (Fs) and (S) are vacuum expectation values of the F- 
component and scalar component of the chiral superfield S that couples directly to the messenger 
sector. Using standard techniques 0] one can then compute the superpartner spectrum and Higgs 
boson mass spectrum. Corrections to the lightest Higgs boson mass are obtained using the full 
one- loop effective potential approximation, as in [12j. (We have checked the MSSM contributions 
against FeynHiggs (26| and we find agreement to within expected uncertainties of 1-2 GeV.) One- 
loop corrections to the pole masses of all strongly interacting particles are also included; these are 
particularly important for the gluino. 

If the exotic states only interacted among themselves and the Higgs fields, then a Z2 quantum 
number could be defined on the superpotential with odd assignments to Q, Q, U, U, E, E and even 
assignments for everything else, leading to stability of the lightest new fermion state. At the 
renormalizable level, the only way the lightest new quark and the r' can decay is by breaking 
this Z2 symmetry via superpotential mixing interactions with MSSM states, 

W mix = e u H u qzU + ^ v H u Qu z - e D H d Qd 3 - e E H d e 3 E, (2.6) 

where e v , e' v , e D , and e E are new Yukawa couplings. Note that this is consistent with matter par- 
ity provided that the supermultiplets Q, Q, U, U, E, E are assigned odd matter parity, so that the 
new fermions have even .R-parity. We assume that the mixing Yukawa couplings are confined to 
the third-family MSSM fields q%, U3,d 3 , £3, in order to avoid dangerous flavor violating effects; the 
bounds on third-family mixings with new heavy states are much less stringent than for first and 



second-family quarks and leptons 27J, |28(. As we will see in section ITVl couplings less than 0.1 to 
third generation quarks and leptons are easily small enough to avoid all flavor constraints. Assum- 
ing this for simplicity, then e v , e' v , e D , and e E are small enough to be neglected in wave function 
renormalizations, and so do not contribute to other couplings' renormalization group equations, 
and only contribute linearly to their own. Furthermore, their effects on the mass eigenstates of the 
new particles can be treated as small perturbations. 

It is interesting to consider the case of S'C/(5)-symmetric interactions near the unification scale. 
If one assigns H u and H d to the 5 and 5 representations respectively, and Q, U, E to the 10 and 
Q, U, E to the 10, then one has 

M Q = M v = M E , (2.7) 
e v = e'jj, e D = e E . (2.8) 

at the unification scale. The further unification in 5*0(10) implies the stronger condition 

e v = e v = e D = e E . (2.9) 

A logical guess is that the origin of the masses Mq,Mu,Me is similar to that of the MSSM \x 
term, and might occur well below the unification scale. For example, one can imagine that they 
arise from non-renormalizable superpotential operators like 

W = -^-SS(X^H u H d + X Q QQ + XuUU + X E EE), (2.10) 
where S, S are SM singlet fields (possibly the same) which carry a Peccei-Quinn charge and get 
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vacuum expectation values (VEVs) at an intermediate scale, as recently proposed in this context 



by [2l|], giving rise to masses \i = X^(SS)/Mp and Mqjj^e = ^q,u,e(SS) /Mp. Note that if the 
dimensionless couplings Xqu,E are small, then their renormalization group evolution from the 
apparent unification scale down to the scale at which S, S get VEVs is the same as that of the 
corresponding masses Mq^u,e, depending only on the wavefunction renormalization anomalous 
dimensions of the chiral superfields Q, Q, U, U, E, E. In this case, it is sensible to evolve the masses 
as if they were the same at the scale of apparent gauge coupling unification, based on an assumed 
unification of the corresponding superpotential couplings \q,u,e- Of course, the relations (|2.7p and 
(|2,9|) are certainly not mandatory. The tree-level relations between couplings (or masses) implied 
by GUT groups can be greatly modified by non-renormalizable terms, alternative assignments of 
the Higgs fields, and mixing effects near the GUT scale. However, eqs. (|2.7|) and (j2.9[) do constitute 
a plausible and useful benchmark case that we will use for some of the explorations in this paper. 
At the TeV scale, typical values obtained from the renormalization group running are then: 

M Q : M v : M E rj 1.8 : 1 : 0.45, (2.11) 

with some variation at the <20% level due to the choice of GMSB messenger scale and A. (The 
ratios Mq/Mjj and My /Me at the TeV scale tend to decrease with larger M mess and A.) The ratios 
of mixing couplings also exhibit a pattern when the unification condition eq. (|2.9p is assumed, but 
with a strong dependence on the trajectory for k. In general one finds e' a slightly larger than e^, 
and e D larger than e v by a factor of 3.5 to 6. In the following, we will sometimes consider the 
typical case 

e v : e' v : e D « 1 : 1.15 : 4.5 (2.12) 

as a benchmark for illustration when considering the branching ratios of tj and b' . 

The model we study here is not the unique extension of GMSB models to include vector-like 
quarks that raise the Higgs mass. One can replace the U + U fields by D + D + E + E fields 
without changing the prospects for perturbative gauge coupling unification, as discussed in [121 ]. 



In that case, a Yukawa coupling H U QD will raise the Higgs mass, and the gross features of the 
superpartner mass spectrum will be unchanged. The exotic fermions will consist of 6' 12 , f, and 



r{ 2 , with decays discussed in 12]. This model is arguably somewhat less motivated, in that it does 
not have complete GUT multiplets. Another variation replaces the 10 + 10 at the TeV scale by a 
5 + 5 + 1 + 1 = L + L + D + D + N + N, with a Yukawa coupling H U LN doing the work of raising the 
Higgs mass. This model has a larger set of possibilities for the GMSB messenger fields consistent 
with gauge coupling unification. However, it also results in a much smaller contribution to Mh, 
unless one includes a larger hierarchy between the exotic leptons and their scalar superpartners. 
In order to keep the present paper bounded, we will not pursue those approaches further here. 



B. Mass spectra for sample models 

In Figure 12.11 we show the mass spectrum of all new particles in a sample model with M mcss = 
1500 TeV, A = 150 TeV, tan /3 = 15, fj, > 0. The left panel shows the result for the minimal 
GMSB model with these parameters, and the right panel the model of interest extended by the 
10 + 10 fields. The minimal GMSB model in the left panel can only manage Mh = 115 GeV, and is 
therefore clearly ruled out if Mh ~125 GeV. In the right panel, we choose Mq = Mjj = Me = 215 
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FIG. 2.1: The mass spectra for new particles in a minimal GMSB model (left) with light messengers 
and a similar model with additional chiral supermultiplets in a 10 + 10 of SU(5) (right). In both cases, 
M mess = 1500 TeV, A = 150 TeV, tan ,3 = 15, and \l > 0. In the model with extra vector-like fields, k = 1, 
and Mq = Mu = Me = 215 GeV at the unification scale. The lightest Higgs mass is 114.9 GeV (left) and 
124.9 GeV (right). 



GeV at the unification scale, as this leads to = 125 GeV. Comparing the two models, we 
see that in both cases M g is close to 1160 GeV; this is significant because the gluino mass is the 
most important parameter pertaining to the discovery of the odd i?-parity sector at the LHC when 
squarks are much heavier, as here. However, in this model at least, the lightest new strongly 
interacting particle is actually the with mass near 700 GeV; it is much lighter than the other 
vector-like quarks b' and t' 2 , and their superpartners, as well as the MSSM squarks and gluino. The 
lightest new particle from the 10+10 sector is the r', which if quasi-stable could also be a candidate 
for the first beyond-the-SM discovery despite lacking strong interactions, as we will discuss below. 
The model with vector-like supermultiplets also produces squarks that are significantly heavier 
than the prediction for minimal GMSB. The Higgsino-like neutralinos and charginos N%, N4,C2 
are also more than a factor of 2 heavier than the prediction of minimal GMSB, corresponding to 
a much larger \fi\. If is treated as a proxy for the amount of fine tuning in the model, we are 
forced to accept that the model with extra vector-like supermultiplets is more unnatural than the 
minimal GMSB model, but this psychological price must be paid if Mh ~125 GeV. 

Figure 12.21 shows a similar comparison, but for a much higher messenger scale M mess = 10 
GeV. The effect of raising the messenger scale is to further increase the squark and slepton masses 
for the model with extra vector-like matter, both in an absolute sense and compared to the minimal 
GMSB model. The Higgsino-like neutralinos and charginos are also much heavier in the extended 
model, pointing to more fine tuning needed in the electroweak symmetry breaking potential, as 
noted above. For the same input parameters, the gluino mass is suppressed in the extended model 
on the right compared to the minimal model, but only by about 4%. In both Figures 12.11 and 
12. 2\ the heavier Higgs bosons ^4°, H°, and have their masses substantially increased when the 
model is extended to include vector-like supermultiplets. 

If [1M2 is positive, there will be a positive correction to the anomalous magnetic moment of the 
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FIG. 2.2: As in Figure [2~2l but with very heavy messengers of supersymmetry breaking at M mcss = 10 14 
GeV. The other parameters are A = 160 TeV, tan/3 = 15, fi > 0, and k = 1 and Mq = M v = Me — 400 
GeV at the unification scale in the model with extra vector-like particles. The lightest Higgs masses are 
115.5 GeV (left) and 124.6 GeV (right). 
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muon, bringing the theoretical prediction into bette r ag reement with the experimental result [29], 
as has been emphasized in the present context by 15j. However, because we are not willing to 
interpret the present ~ 3<r discrepancy as evidence against the SM, we simply take 11M2 > and 
do not impose any constraint from (g — 2)^. It is also useful to note that for all models of this type, 
the effect of the vectorlike quarks is to bring slightly closer agreement with precision electroweak 
oblique corrections than in the SM, but not by a statistically significant amount [l2| . 



C. Achieving M h -125 GeV 

The corrections to the lightest Higgs mass are most strongly dependent on the masses of t^t^ 
and their superpartners ^234; with AM^ increasing with the hierarchy between the average scalar 
and fermion masses. The masses of 2 3 4 scale with the super symmetry-breaking parameter A, and 
the smaller they are, the smaller the fermion masses t'i 2 and b' must be in order to accommodate 
M/j ~125 GeV. The masses of the gluino and t[ are of particular interest, since pair production of 
one of them is likely to give the initial discovery signal at the LHC. Figure [231 shows (green sloped 
funnel) regions in the M t ^ vs. Mg plane in which 122 GeV< Mh <128 GeV, for tan/? = 15, with 
k = 1 at the unification scale. The variation in M g is obtained by varying A, and that of Mjv by 
varying Mq = M\j = Me at the unification scale. Three choices of the messenger scale are shown, 
M mess = 10A, 10 10 GeV, and 10 14 GeV. Note that, pending exclusions by direct searches for gluino 
and ij, it is easy to obtain M^ ~ 125 GeV in this class of models, with M g lower than 700 GeV 
and M t i lower than 300 GeV even if the messengers are light. Therefore, each new search result 
at LHC probes an interesting region of parameter space consistent with M^ ~125GeV, unlike in 
the usual GMSB models. 

The dependence on tan f3 is shown in Figure [231 with allowed regions for 122 < Mh < 128 GeV 
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FIG. 2.3: Regions in the M t > vs. Mg plane for 
extended GMSB models with a minimal mes- 
senger sector and a 10 + 10 of vector-like quarks 
and leptons, with M mcss = 10A and 10 10 GeV 
and 10 14 GeV. The vector-like fermion mass pa- 
rameters Mq and Mjj are taken to be unified 
at the gauge coupling unification scale. The 
Yukawa coupling k is near its fixed point, ob- 
tained by setting k = 1 at the unification scale, 
and tan/3 = 15. The green sloped funnel re- 
gions have a calculated Higgs mass satisfying 
122 GeV < M h < 128 GeV. In the orange 
shaded region, M T > < 100 GeV, assuming that 
Me = Mq = Mjj at the unification scale. 



in the tan /3 and M g plane. In each graph, the lighter green curved region furthest right corresponds 
to the choice of M t i = 1200 GeV and the darker green curved region to the left of it corresponds 
to Mtf = 600 GeV. The upper left triangular red region corresponds to M T i < 100 GeV. The three 
graphs shown correspond to M mess = 10A and 10 10 GeV and 10 14 GeV, and all have k = 1 at the 
unification scale. More details regarding underlying parameters are found in the figure caption. 

Note that an intermediate value of 10 < tan/3 < 35 enables a lighter gluino mass, and so 
lighter MSSM squark masses, than found for tan/3 outside of that range. For larger tan/3, the 
corrections to from the tau-stau sector are negative and big,^ so that larger supersymmetry 
breaking masses (indicated in the plot by Mg) are required. For tan /3 < 10, the tree-level Mh is 
much smaller, requiring heavier superpartners to obtain M^ ~125 GeV. Similar figures are found 
in [rol . [2ol | . but with Mq = Mjj at the TeV scale, rather than at the unification scale as chosen 
here. An important point 20j is that there is an upper bound on tan/3 in these models, following 



The tau-stau loop contributions are larger than the bottom-sbottom ones, despite having a smaller Yukawa coupling 
and no color factor, because the staus are much lighter than the sbottoms. 



9 



M = 10 A 

mess 



M ™« = 10 10 GeV 

mess 




1000 1500 2000 
Gluino mass (GeV) 

J4 



_i I I l 



1000 1500 2000 
Gluino mass (GeV) 



2500 





60 - 




50 1 




40? 


ca 




c 

03 


30 : 


■*—> 






20: 




10 L 



M = 10' 

mess 



GeV 




FIG. 2.4: Regions in the tan/3 vs. M g plane 



14 



00 1000 1500 2000 
Gluino mass (GeV) 



j i i i_ 



2500 



with M mcss = 10A and 10 lu GeV and 10 
GeV. The vector-like fermion mass parameters 
Mq,Mu are taken to be unified at the gauge 
coupling unification scale. The Yukawa cou- 
pling k is near its fixed point, obtained by set- 
ting k = 1 at the unification scale. The green 
shaded curved regions have a calculated Higgs 
mass satisfying 122 GeV < M h < 128 GeV, 
with the darker green region corresponding to 
Mf = 600 GeV and the lighter green region 
having M t / = 1200 GeV. In the red shaded re- 
gion, M? t < 100 GeV. Below (above) the dashed 
line the NLSP is a neutralino (stau) . Above the 
solid curves, the vacuum is unstable by tunnel- 
ing to a state with VEVs for the staus, with a 
liftime shorter than the age of the universe. 



from the general bound obtained in [30)] by requiring the standard electroweak-breaking vacuum 
to be stable (with a lifetime longer than the age of the universe) against tunneling to a vacuum in 
which the stau fields have VEVs. We show this bound for our models as the solid lines in Figure 
12.41 We see again here in this figure that gluino masses easily accessible by LHC now or in the near 
future are sufficient to deliver a light Higgs boson of mass ~125 GeV, and this can be achieved for 
Mg < 2.5 TeV even if tan f3 is as low as about 3. 

As remarked above, Mq and M\j are independent in a general theory. Figure 12.51 explores 
this freedom by showing lines in the My , My plane that predict Mh = 125 GeV, for models with 
tan/3 = 15 and M mess = 10A, with the ratio Mq/Mjj allowed to vary. The special cases with 
Mq = Mjj at the unification scale (as in the examples of Figures [2.11 12. 2\ I2.3[ 12.41 above) are noted 
by green stars. The three curves correspond to A = 120, 180, and 240 TeV, resulting in M g ~ 960, 
1390, and 1800 respectively. (There is some small variation in the gluino masses on each curve.) 
We find that for equal values of other parameters, Mh remains approximately constant for fixed 
values of the arithmetic mean of M# and M t i . In particular, the geometric mean is not as good 
a figure of merit. For each curve in Figure 12.51 we see that there is no minimum value of M t i 
from the Mh ~125 GeV constraint alone, because one can always take a very large or small ratio 



10 



- 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 — 

M~ = 1800 GeV 




M h = 125 GeV 



IV1 = 10A 

mess 



500 1000 1500 2000 
M t , (GeV) 



FIG. 2.5: Lines in the M t ^,M t i plane with pre- 
dicted Mh = 125 GeV, for models with tan j3 — 
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of Mq/M{j. However, on each curve corresponding to a fixed Mg, the requirement ~125GeV 
implies a minimum value of M t i , and a maximum value of M t ' . 



D. Comment on gravitino dark matter 



In GMSB models, the LSP is likely to be the gravitino G, with mass Mq = AM mcss /\/^Mp, 
where M p = 2.44 x 10 18 GeV. In principle, the gravitino could be a dark matter candidate. One 



possibility is the gravitino superwimp scenario 3l| in which the gravitino abundance is assumed to 
be suppressed by a low reheating temperature or diluted by some other non-standard cosmology, 
followed by the bino-like neutralino LSP freezing out and then decaying out of equilibrium according 
to Ni — > jG, with a lifetime given approximately by [32] 



Tjj = 7.5 x 10 sec 



G 



GeV 



100 GeV 



M 



(2.13) 



If N\ — > ZG is kinematically allowed, then this lifetime is reduced by a factor 1 + 0.3(1 — m^/m^ ) 4 
[33)]. If the gravitino is to be a significant component of the dark matter, this lifetime should be 
smaller than about 0.1 to 1 sec, in order that the successful predictions of primordial nucleosynthesis 
are not affected. This is in tension with a cosmologically relevant relic abundance of gravitinos 
from decays of thermal binos, given by 



n& h 2 = ^n^h 2 . 



m 



Ni 



(2.14) 



Here 



n^h 2 ps 0.013 ^ ^ 1 



r(l + r 2 ) VlOO GeV 



[l + 0.07 ln(v^ 100 GeV/mgJ] 



(2.15) 
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FIG. 2.6 

1 second, 0.1 second, 1 meter, and 1 cm (solid 
lines), and gravitino abundance resulting from 
decays of thermal neutralino NLSPs tt^h 2 — 
0.11, 0.01, and 0.001 (dashed lines), in the rrig 
and M mess plane. The gluino mass variation 
was obtained by varying A, with tan/3 = 15 
and k = 1. 



is the relic density of binos that would be found today if they were stable, given in a convenient 
approximation 3J|, with r = m 2 ~ /rn~ . To illustrate this, we show in Figure [2U1 solid lines of 



constant = 0.1 and 1 second (relevant for nucleosynthesis) and 1 cm and 1 meter (relevant 



for collider physics), compared to dashed lines of constant Q^h 2 = 0.11, 0.01, and 0.001, in the 
(rrig, M mess ) plane, with the variation in gluino mass obtained by varying A. It is difficult to rec- 
oncile gravitino dark matter with the standard picture of primordial nucleosynthesis in this model, 
without going to very large superpartner masses (rrig 3> 2.5 TeV), in which case the vector-like 
quarks would not be necessary and prospects for any discovery of new particles beyond h° at the 
LHC would be exceedingly grim. J Such a massive superpartner spectrum runs counter to the pur- 
pose of this paper, which aims to accommodate the ~ 125 GeV with lighter superpartners accessible 
to the LHC. In the scenario considered in the present paper, these considerations suggest that dark 
matter is composed mostly of axions or some other particles, with a negligible contribution from 
gravitinos, and messenger mass scales much above 10 11 GeV are therefore apparently disfavored 
as indicated in Figure [ 



III. MASSES OF EXOTIC QUARKS AND t[ DECAYS 



Taking into account the full superpotential of the theory Wmssm + Wqjje + Wmix the fermionic 
mass matrices for up-type and down-type quarks are [H] 



'VtVu tijVu ^ 

Mu k'v d 
\^ v v u kv u MqJ 



Mr 




Mr 



(3.1) 



with mass eigenstates t, t'^t'^ and b,b' respectively. The zeros appear as a consequence of a choice 
of basis. As mentioned earlier, we assume that e u , e^, and e D can be treated as small perturbations 



* Two recent papers [35l . |36| ] have noted the complementary approach that in normal gauge mediation models, one 
can accommodate gravitino dark matter and Mh ~125 GeV, at the cost of such very heavy superpartners. 
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FIG. 3.1: Masses (left panel) and t[ branching ratios (right panel) for vector-like quarks. Here, M? is fixed 
to 600 GeV, and the ratio of mass parameters Mq and Mjj at the TeV scale is varied. For small (large) 
Mq/Mjj, the state t[ is mostly SU(2)l doublet (singlet). The green band shows the region obtained with 
the unification condition Mq = Mjj imposed at the gauge coupling unification scale. The left edge of this 
band corresponds to M mess = 10 14 GeV, and the right edge to M mcss = 10A. The weak-scale parameters e v , 
e' n , and e D that describe mixing of the vector- like quarks with the top and bottom quarks are in the ratio 
1 : 1.15 : 4.5, which are typical approximate values predicted by requiring them to be unified at the gauge 
coupling unification scale. 

in these mass matrices. Then one always finds M t t < My < M^, and the exotic quarks will decay 
according to t' 2 -> Wb', ht[, Zt' x and b' -> W^h[, Wt, Zb, hb and t[ -> Wb, ht, Zt, when kine- 
matically allowed. Formulas for these decays widths, which will be used in our phenomenological 
discussion below, can be found in Appendix B of [12j, and in a more general framework in the 
Appendix of the present paper. 

In Figure 13.11 we plot the mass eigenvalues of the exotic quark states t'i 2 and b' as a function of 
Mq/Mjj in the left panel, and the branching fractions of t^ vs. Mq/Mjj in the right panel. Within 
this figure m t i is fixed to be 600 GeV. For Mjj <C Mq the t[ state is a nearly pure ££7(2)£-singlet, 
and it decays into Wb, ht and Zt primarily through the interaction e u H u q-^U. The dominant decay 
mode in that limit is to Wb at slightly over 50%, but ht and Zt final states are non- negligible. In 
the opposite limit Mq <C Mjj, the state t\ is nearly pure 5J7(2)^-doublet, and it decays mostly 
into ht, with Zt a significant subdominant mode. Note that the case Mq ~ Mjj at the TeV scale 
is actually in a transition region for the branching ratios. These results were obtained assuming 
that e v , e'jj and e D are in the low-scale ratios of 1 : 1.15 : 4.5, which are approximate results 
from assuming they are unified at the gauge coupling unification scale. The thick vertical band in 
Figure [3~T1 indicates the ratio of Mq/Mjj at the TeV scale under the assumption that Mq = Mjj 
at the gauge coupling unification scale (typically in the range ~ 3-8 x 10 16 GeV for these models). 
The left edge of this band corresponds to M mess = 10 14 GeV, while the right edge of the band 
corresponds to M mess = 10 A. 

It is also interesting to consider the dependence on the mixing couplings t u , e' v , and e D , because 
the relation eq. (|2.12p may not hold. This is illustrated in Figure 13. 2| in which we hold fixed 
Mq/Mjj = 1.8, and vary e'^/e^ with e D = 0, and e D /e v with e' a = 0. When the ratio /e^ | is less 
than a few, and when e D /e v < 50, one recovers results similar to the unification-motivated results 
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FIG. 3.2: Branching ratios for t[ as a function 
of ratios of the mixing couplings. In the top 
two panels, the ratio e^/e^ is varied over posi- 
tive (left panel) and negative (right panel) val- 
ues, with e D — 0. In the lower panel, the ratio 
e D /eu is varied, with e' v = 0. (The results do 
not depend on the sign of the ratio in the latter 
case.) 
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given in Figure I3.ll This is because in that case the effects of e v are dominant because of the 
S'{7(2)L-singlet nature of t[. However, for larger values of e'^, one enters a 'W-phobic" regime for 
ti in which the ht final state can dominate with B(Wb) very small. Conversely, for e D very large, 
one goes over into the charged-current dominated case that B(VF6) = 1, which coincides with the 
prediction for a sequential t', the subject of most experimental searches. Clearly, it is crucial that 
experimental searches go beyond this case, to take into account and hopefully exploit the ht and 
Zt final states. 

The dependence of these branching ratios on the magnitude of the t\ mass is mild provided that 
it is well above the W, Z, h masses. This is illustrated in Figure 13.31 which shows the branching 
fractions of t^ as a function of its mass, keeping fixed k = 1 and using the unified boundary 
conditions Mq = 1.8Mu and e L , : e' v : e D = 1 : 1.15 : 4.5. For low t[ masses < 400 GeV, the 
branching fractions show some variation, but with higher t^ mass they asymptote to B(VF6) = 50%, 
B(ht) = 25%, and B(Zt) = 25%, but with B(ht) > B(Zt) for finite masses relevant to the LHC. 
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FIG. 3.3: Branching ratios for t[ as a function 
of its mass, obtained with Mq — 1.8Mv and 
e v : e' v : e D = 1 : 1.15 : 4.5. 



IV. PRECISION TESTS FROM MIXING WITH THIRD-FAMILY FERMIONS 

The introduction of an additional b' quark that mixes with the third generation b quark can 
induce a tree-level shift in the Z boson coupling to the right-handed b quark mass eigenstate 
compared to the SM. Such a shift is very severely constrained by the measurement of R b at 
LEP [37], with 

R ex P t = o 21629 ± 0.00066 (4.2) 
from {38|. The SM computed best fit value is 3^, 35^] 

R^ M = 0.21579 ± 0.00013. (4.3) 
Thus, the 3a range of allowed shifts in Rb compared to the SM value is 

- 0.0015 < 5R b < 0.0025, (4.4) 
where 5R b = Rh - R? M . From eqs. (|A.9p and (jA.lOp in the Appendix, and relating the coupling 
tree-level shifts in the couplings are 



conventions in [38|, [39J to ours by g L = (c-w 1 9)9 ,\ , and g R = —{cwlg)g^- , we see that the 

^ Q 3 3 



- i> , h l lr ,/ 12 e 2 D v 2 cos 2 /3 

6g b L = 0, 5g b R = --R' 4 3 2 &- ^ , (4.5) 

I Imp 

which shows that the mixing always reduces the magnitude of the right-handed b quark couplings 



to the Z boson. With this definition the resulting shift in R b is 37] 



R b = Rl M (l + 0M5 5g b R ) (4.6) 
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which implies that the 3a range of allowed 5g b R is 

- 0.011 < 5g b R < 0.018. (4.7) 

Thus the requirement that Rb is in 3-0" agreement with experiment gives a constraint on e D , tan/3 
and my. From eqs. (|4.5p and (|4.7p we find the requirement that 



leJ < 0.42 tan/3 ( Wi + — L— . (4.8) 

! \500GeV/ V tan 2 /3 v ' 



A similar analysis follows from considering shifts in Ab and A b FB . The current 38|, [39] experi- 
mental situation is that 

A expt = o g23 ± Q 02Q A b£cpt = Q 0Q92 ± q q 016) (4 9) 

whereas the SM computed best fit values are 

Af M = 0.9346 ± 0.0001 and = 0.1033 ± 0.0008. (4.10) 

Let us focus on A b FB , as the SM prediction is 2.8a too high compared to the measurement (see 
table 8.4 of 13). 

From the definition Afb = (#£ — 9r)/(9l + 9r) one can compute the shift in A b FB from a shift 
in g b R to be 

4 B = ^ S /(1-1.74) M^ = -0.18^. (4.11) 

Since < 0, this implies that the shift in the prediction of A b FB is always positive, increasing 
the tension between theory and experiment. If we therefore assume that the b — b' mixing is no 
more than a la effect in the "wrong" direction (i.e., SA b FB < 0.0016 from b — b' mixing), this puts 
a limit on 8g b R that translates to exactly the same formula as eq. fj4.8j) except that 0.42 is replaced 
by 0.38. Thus, the constraints on b' mixing are not very severe as long as my is greater than a few 
hundred GeV or tan f3 is not small. 

Another way to constrain the mixing of SM third-family quarks with the exotic quarks is through 
the CKM matrix element Vtb- Here, we cannot assume unitarity of the CKM matrix, since it will 
not be in general [see eq. (|A.13|) in the Appendix]. If the e D coupling is present simultaneously 
with the e v or couplings, then the situation is complicated by the fact that the W boson will 
have small couplings to right-handed SM quarks as well as left-handed quarks. For the sake of 
illustration, consider the case that only e u is important, and suppose that the SM Yukawa coupling 
matrices are such that if e v were exactly 0, then Vtb would be very close to 1 (as one finds in the 
SM with CKM unitarity assumed), so that all mixing of the first two families with the third family 
and the vector-like quarks can be neglected. With those assumptions, from eq. ()A.13P we obtain 



500 GeV\ 



1-^^0.06 4si^/3( ^ \ , (4.12) 
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This can be compared to the values obtained from single top production without assuming CKM 
unitarity, V tb = 0.88 ±0.07 (from Tevatron 0) and V tb = 1.04 ±0.09 (from CMS [41]). Thus, even 
if e v is near unity and t\ is not much heavier than its experimental bound, the CKM constraint 
does not impact the model. 

Next, we consider the implications of r — t' mixing. This mixing will induce a positive shift 
in the g T L = {g/cw)g Z ^ coupling to the Z boson, while g T R = — {g/cw)g%. is unaffected. Prom 

eqs. (|A39|) and (TOoj . 



^ = ^ |2si 4 lf ^ = °- (413) 



An important effect that results from this shift is an alteration in the A T observable. From the 
definition Ai = (g\ — g\)/{g\ + g\), the shift in 5A T from a shift in 5g T L is 



A r = A* M (1- 23 5gl) (4.14) 

which demonstrates the high sensitivity to changes in the r lepton couplings to the Z. 
The experimental and theoretical values [38] of A T are 

A e T xpt (P T ) = 0.1465 ± 0.0033 and A S T M \P T ) = 0.1480 ± 0.0011. (4.15) 

Keeping the prediction to within 3a of the experimental measurement requires that —0.0120 < 
5A T < 0.0090. Since 5A T is always negative from the r' mixing, the lower limit is the applicable 
constraint. From eq. (|4.14p we see that 5g T L < 0.0033, or 



l £ -l <0 ' 23tM ^(50?§ev)V 1 + i^' ( " 6) 

This requirement is not terribly constraining, especially considering that the SM r Yukawa coupling 
y T = 0.01 is much smaller than the general constraints on €e when m T i > 100 GeV. 

Finally, one can attempt to constrain the t-t' mixing through the r decay measurement. The 
analysis of [13] corresponds to | C/43 1 2 < 0.0053 in the notation of the Appendix of the present paper, 
which therefore implies the same constraint as eq. (|4,16p but with 0.21 replacing 0.23. However, 
this is a 1-a constraint. Also, this assumes that the PMNS matrix is unitary, and that mixing in 



the electron and muon sectors is absent, which need not hold 43J]. In any case, there is no impact 



on the coupling e E in this model unless tan j3 is small, and the r' is light. 



LHC PHENOMENOLOGY 



The exotic quarks could in principle have a significant effect on the production and decay of 
the lightest Higgs boson. For an additional chiral fourth family, which relies entirely on Yukawa 
couplings for its large masses, there is a very largepositive effect on the production cross-section 
0, 44 1, in strong conflict with the current limits [H, 0]. However, in the vector-like model under 



present consideration, the situation is very different. The corrections to the hgg and /177 effective 



interactions can be found from the general formulas in 45[]. Applying these, we find that for the 
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Gluino Mass (GeV) Gluino Mass (GeV) 

FIG. 5.1: NLO production cross-sections for gg and gQ + gQ and CfN\ and C\G\, as a function of the 
gluino mass, for GMSB models with extra vector-like quarks in this paper with M mcss = 10A and tan/3 = 15, 
in pp collisions with ^fs — 8 and 13 TeV. 



case k! <C k ~ 1, these corrections are totally negligible. Even if k' is sizable, the corrections to 
gg — > h — > 77 are quite modest, at most at the 5% level for M t i = 500 GeV and tan f3 = 5 and 
k = 1, k' = 0.7, and can have either sign depending on the relative phases in the t^t^b' sector 
mass matrices. For larger tan /?, the size of the effect decreases. We conclude that at least for LHC 
physics in the short term, the loop effects of the exotic quarks on Higgs production and decay are 
probably too small to hope to observe. 

The model under consideration differs from other variants of the MSSM in that there are two 
distinct paths to a new physics discovery. First, we may discover the odd i?-parity superpartners 
of the SM states. Second, we have the exotic quark and lepton states. These two possibilities are 
essentially decoupled, and it is unclear which of them should provide the initial discovery of physics 
beyond the SM, since the masses and decays are negotiable within the general model framework. 
We will begin by commenting on features of the superpartner phenomenology at LHC, making the 
comparison to other standard searches. 



A. Superpartner signals 

If the NLSP is a neutralino (N\) that is stable on detector-crossing time scales, the resulting 
phenomenology is very similar to "standard supersymmetry" signatures (e.g., mSUGRA). The 
squarks are comparatively heavy, with up and down squarks, which play the most important role 
in LHC production, between about 1.6 and 2.3 times heavier than the gluino (see for example 
Figures 12.11 12. 2\i . Therefore, the discovery potential comes mostly from gluino pair production, 
gluino+squark production, or the production of wino-like charginos and neutralinos, followed in 
each case by decays to jets, leptons and large missing energy. The production cross-sections 



computed to next-to-leading order by Prospino 46] are shown in Figure [57T1 for the most important 

processes pp — > gg and gQ + gQ and C^iV^ and Cf Here we used a model line with M mess = 
10A, tan/3 = 15, k = 1, and \x > 0, but the dependence on these particular assumptions is 
mild, with the exception of gQ + gQ, which becomes smaller for a given M g if M mess is larger. 
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Although the gluino+gluino and gluino+squark pair production cross-section are smaller than the 
chargino-neutralino rates for M g > 650 GeV at y/s = 8 TeV, and for M g > 1050 GeV at yfs = 13 
TeV, the gluino and squark signals should have higher acceptances due to more visible energy. 
However, any attempts to probe much beyond M g = 1 TeV at y/s = 8 TeV may have to rely on 
chargino/neutralino production rather than gluino/squark production. 

The branching ratios of the gluino are, for the typical low-M mess model in Figure 12.11 

g -> jjCx (38%), tbCt (17%), jjN 2 (19%), bbN 2 (6%), 

ttN 2 (3%), jjNi (12%), bbNi (1%), ttN x (4%), (5.1) 



from SDECAY 47(, where j denotes a jet from a u,d,s,c quark, and the notation omits the 
distinction between quarks and antiquarks. Up and down squarks essentially always decay to a 
gluino and a very energetic jet. The wino-like charginos and neutralinos decay almost entirely 
through the lighter stau, which then decays as f± — > tN\ with a branching ratio of 100%: 

Ci -> f x v^TvN x (96%), WN X (4%), (5.2) 
N 2 -> nr^T+r-Ni (96%), hNi (4%), (5.3) 



Thus a high proportion of events will have 2, 3, or 4 taus in the final state, manifested either as 
hadronic tau jets or softer e, fi. This is an important difference compared to mSUGRA, where 
comparable models with such heavy squarks have large rriQ and therefore also have heavy staus, 
and so cannot produce such a predominance of taus in the final state. 

In contrast, models with higher M mess will have > Mg i « , as illustrated by the 
example in Figure 12.21 implying a much lower tau multiplicity. In that example model, we have 
for the gluino decays 

g -> jjCx (33%), tbCi (18%), jjN 2 (16%), bbN 2 (6%), 

ttN 2 (3%), jjNi (13%), bbNi (2%), ttNi (9%), (5.4) 



similar to eq. (|5. 1 1) , with a slightly higher average number of b jets. However, the wino-like charginos 
and neutralinos decay very differently than in the low-M mess case: 



C x -»• WN X (100%), (5.5) 
N 2 -»• hN x (97%), ZN X (3%). (5.6) 

This means that over 40% of gluino pair production events, and almost all C\N 2 events, will have a 
Higgs boson in them. For such models with Mf ± > ~ > the signals are sufficiently similar 
to mSUGRA ones with large mo that one can safely approximate the limits by those obtained by 
ATLAS and CMS for the same gluino mass and heavier squarks. The ratios of squark masses to the 
gluino mass are shown for our model in Figure [5T2l These squark/gluino mass ratios correspond 
approximately to CMSSM models with mo/Mi/ 2 ranging from about 3.4 (for low M mess ) to 5.2 
(for high Mp^). At present, the LHC limits for these large tuq cases imply only M g > 850 GeV 



from [48l. |49| . A direct comparison is hindered somewhat by the fact that the LHC collaborations 
unfortunately choose to present results for the CMSSM in terms of the unphysical input variables 
(rriQ, Mi/ 2 ) rather than physical gluino and squark masses. 
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FIG. 5.2: The ratio of the first-family 
squark masses to the mass of the gluino, 
as a function of M mcaa . Here, tan/3 = 15, 
H > 0, and A = 160 TeV. 



FIG. 5.3: The ratio M f jM 6i of the lighter 
stau and chargino masses, as a function 
of M mess , for tan/3 = 10,20,30,35,40,45. 
Here, A = 160 TeV and fi > 0. The missing 
parts of the lines for tan/3 = 35,40,45 are 
due to the vacuum stability requirement. 



Because of the importance of the transition in parameter space between the cases that f\ is 
lighter or heavier than the wino-like neutralinos and charginos, we show in Figure [531 how the ratio 
M? 1 /Mg behaves as a function of M mess for various values of tan /?. For Mf^jMg < 1, the decays 
C\ — > f\V and Ni — > t\t dominate; otherwise, decays to W and h dominate. Depending on tan/3, 
we see from Figure 15.31 that the transition between these two regimes occurs at an intermediate 
scale of a few times 10 9 to a few times 10 11 GeV. 

If the decay N\ — > 7G is prompt, then the above event topologies will be supplemented by 
two energetic isolated photons, for which SM backgrounds are quite low. This would increase the 
discovery potential dramatically, and would probably guarantee that the discovery would happen 
in the C1N2 production channel, due to its larger cross-section. Because the NLSP decay width 
is proportional to 1/F 2 , where F > AM mess , we see from Figure [5T3l that the prompt neutralino 
NLSP decay signal should be tttjj + E™ ms , where r can be either a softer lepton or a hadronic 
tau jet. 

Another possibility is that the NLSP is the lighter stau, which can only occur in our model 
framework if tan f3 is large. (However, tan (3 cannot be too large, and M mess must be low, given the 
constraints on vacuum stability evident in Figure [2~4"1 ) In that case, all superpartner decays chains 
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will terminate in f\ —> tG, where G is the goldstino (gravitino). In each decay chain from a gluino, 
chargino, or neutralino parent, lepton flavor conservation dictates that there is another r produced. 
This means that if the NLSP stau decay is prompt, essentially all super symmetric events will have 
at least 4 taus, while if it is not prompt, one has at least 2 taus and 2 quasi-stable staus which 
can be detected as slow-moving heavy charged particles. However, the parameter space in which 
there is a stau NLSP is limited, as one runs into the constraint from stability of the vacuum noted 
2(3, 30]. For example, in the models of Figure l2~4l one sees that this requires a low M mess , and 



m 



Mg < 1100 GeV, and 21 < tan/3 < 34. 



B. Search for t[ at the LHC 

The production cross-sections for generic exotic heavy quarks at the LHC are shown in Figure 
15.41 for various yfs values. The collider phenomenology of the t\ depends crucially on whether 
it decays promptly or not. If the mixing between the exotic quarks and the SM quarks is very 
small, then there is a chance that t' x could be stable on time scales relevant for collider detectors. 
Assuming first the unification ratio Mq/Mjj = 1.8, so that is mostly SU(2)l singlet, we find a 
lifetime for t'i of 

0.94 mm for e = e v sin (3, e' a = e D = 0, 
290 mm for e = e' v sin /3, e a = e D = 0, (5.7) 
340 mm for e = e D cos /3, e v = i v = 0. 

For simplicity, we have taken the limit M$ 3> in eq. (|5.7p . For masses closer to the weak boson 
masses, kinematic factors increase the lifetime somewhat. To illustrate the opposite limit of the tj 
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FIG. 5.5: Limits on m f ; vs. — > bW), from 
the CMS search for the final state bbWW based 
on 4.7" 1 fb reported in [53]. 
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being mostly an SU(2)l doublet, consider Mq = 0.5Mu, which results in 



cr 
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(5.8) 



If we require for the definition of prompt decays that cr < 1 mm, then we need only either e v or 
e' a to be greater than a few times 10~ 7 to ensure prompt decays. Note that the e D contribution to 
the inverse lifetime is suppressed by cos 2 j3. 

Let us first assume the case of decaying promptly. The LHC experiments have several analyses 
based on the production of heavy top-like quarks. The most stringent direct search bounds are 
from CMS, but are limited to the extreme cases that either the Wb or the Zt final state dominates. 
For B(t' -)■ Wb) = 1, CMS obtains M t , > 557 GeV using 4.7 fb" 1 
they obtain M t < > 475 GeV using 1.14 fb" 1 



5l|, and for B(t' -> Zt) = 1, 



52|. However, in our case the branching ratios are 



split among the final states depending on the mixing couplings, as seen in Figures f3.HI3.3i In much 
of parameter space, where M t i > few hundred GeV and Mq > Mjj, we find Bit^ — > bW) ~ 50%. 
Therefore, requiring t'-^-y — > bWbW means a reduction by a factor of 4 in total cross-section 
applicable for the analysis. Taking this factor into account, and comparing the cross-section limits 
at LHC as derived in [53|] with the total direct rate in our theory assuming B(t[ — > bW) = 50%, 
we extrapolate to find the current limit to be M f / > 420, even without using the other final states. 
This is consistent with another recent analysis 5J] . In Figure 15.51 we show the limits as a function 
of B(t± — > bW) based on the t'\i'-y — > bWbW limits only. Of course, a more general search using all 



three final states Wb, Zt and ht will find a stronger limit. In ref. 54 |a reanalysis of these direct 
search limits together with a reinterpretation of an ATLAS search 55J] for b' —> Wt in terms of 
t' pair production is argued to give a bound M t / > 415 GeV, for any combination of the three 
branching ratios for t' — > Wb, Zt, Zh. Going forward, the detector collaborations should strive to 
incorporate all three final states in their search strategies as much as possible, in order to maximize 
the model-independent reach in the t[ mass. For any value of M t i , the mixing couplings can be 
chosen in such a way that any of the Wb, Zt, or ht is the dominant decay mode, and they may 
all be comparable to each other. This should be kept in mind in the planning and interpretation 
of hadron collider searches. Even if t' — > Wb has the largest branching ratio, searches with mixed 
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final states (t' — > Wb){t' — >• Zt) or (t' — > Wb)(t' — > ht) may give the strongest signal, exploiting the 
presence of Z — > £ + £~ and 2, 3, or 4 fo-tagged jets, or even h — > bb or h — > t + t~ with a "known" 
invariant mass of ~125GeV. This is especially important given that there are other, compl etely 
different, new physics models that predict exotic quarks within the reach of the LHC [561]- [641] . 
which can span the possible branching ratios into these three final states. It would be especially 
interesting to observe and study events with h — )■ bb or h — > t + t~ in tj production, since these 
decay modes are quite difficult to observe at the LHC from direct Higgs production. 

If the t\ is stable, it can be searched for as a strongly interacting heavy stable charged particle. 
The implications for the search are expected to be similar to that of a quasi-stable top squark when, 
for example, it is the NLSP and the decay to gravitino is very suppressed and the lifetime is greater 
than the size of the detector, cr > ^detector- The search strategy [65f] relies on first identifying large 
dE/dx energy depositions in the inner tracker due to the massive stable charged particle traversing 
it. This combined with the requirement of high px is the so-called tracker method of discovery. 
In addition the excellent timing of the muon system enables a time-of- flight cut, since a massive 
particle will have smaller velocity usually than a muon and thus takes more time to reach the 
outer muon chambers. The combination of these two methods, tracker and time-of- flight, yields 
powerful constraints from the ^fs = 7TeV data. With 4.7 fb -1 of integrated luminosity, we can 
compare the cross-section vs. mass limits of [53J] to the cross-section computation in Figure 15.41 
and from extrapolation of these results conclude that there is a limit of quasi-stable mass of 
m t' 1 ~ 950 GeV. We estimate that more than twice this sensitivity could be achieved at 14TeV 
LHC with more than 10 fb -1 of integrated luminosity. 



C. Search for b' at the LHC 

In addition to the t' x , the 10 + 10 model has exotic quarks b' and t' 2 - It is of particular interest 
to ask what are the sensitivities to b' production at the LHC [fjf]], since its mass may be nearly 
that of the t\ fermion when Mq < Mjj, as seen in Figure [37T1 Given that its production rate is 
nearly the same as that of a similar mass t^, due to QCD contributions dominating, we must ask 
how the LHC would find this state, which is almost pure 5C/(2)i-doublet in both its right- and 
left-handed components. 

The b' can have two-body decays through the mixing parameters e v , e' a , or e D to possible final 
states Wt, Zb and hb. Again, we are assuming that the exotic fermions couple only to the third 
generation weak eigenstates in order to tame potential flavor problems in the theory. The decay 
widths are calculated (in the more general case of arbitrary mixing with the SM quarks) in the 
Appendix, eqs. ()A.15|) - (jA.17p . The relative fraction of the decays into Wt versus Zb and ht depends 
to a large extent on the ratio e^/^ tan/3). If this ratio is smaller than 1, or if e D /e v is small, 
then b 1 yields mostly Wt, if kinematically accessible. If the ratio is larger than 1, then the b' yields 
mostly Zb and hb. The branching ratios are shown in Figure [5TU1 for a nearly pure doublet b' , as a 
function of e D /eu with e'^ = l.le^ and e v = 0. 

It is also necessary to consider the flavor-preserving decay b' -> W^t^. If M v > M t > + M w , 
then this is an on-shell two-body decay and it will dominate. However, for the case that b' is mostly 
doublet, the decay will be three-body with the W boson off shell. The formula for this decay width 
is found in the Appendix, eq. (|A.18|) . In Figure IBTTl we show this width for the idealized case 
that b' has pure doublet couplings to W and t! x , as a function of the mass difference My — M t > , 
which is the most crucial parameter. For comparison, the two-body flavor-violating decay widths 
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Mjj = 600 GeV, 
M f = 592 GeV, 
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FIG. 5.6: Branching ratios for a nearly pure 
SU(2)l doublet b' into Wt, Zb, and hb, as a 
function of e D /tu. Here Mq — 600 GeV and 
Mu = 1200 GeV with k = 1, so that m 6 ' = 600 
GeV and rrif = 592 GeV. The solid lines have 
e„ = l.le U} and the dashed lines have e' v = 0. 
The three-body decay b 1 — > t[ff through an off- 
shell W boson is highly suppressed by kinemat- 
ics, and is assumed to have a small branching 
ratio. 



FIG. 5.7: The decay width T(b' ->• W^h[) 
as a function of the mass difference My — 
M t > , assuming that b' and t[ form a nearly 
unmixed SU{2) L doublet, for M v = 500 
GeV. (The results for the range 400 GeV 
< My < 1000 GeV are visually nearly in- 
distinguishable from the line shown on this 
graph.) 
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are approximately: 



0.1 GeV 



1000 GeV / V0.1, 



(5.9) 



for the cases e = e' v sin j3, e D = e v = and e = e D cos /3, e v = e'^ = 0, and 



9 x 10~ 5 GeV 



My 



1000 GeV 



1000 GeV J V Mu 



) 4 (^) ; 



(5.10) 



for e' v = e D = 0. Thus, the decay b' W^t^ may or may not dominate in this case, with a strong 
dependence on both the mass difference and the mixing couplings. 

If b' mostly decays into Wt the current limits arise from a search by CMS [67] based on 4.9 fb _1 
of integrated luminosity, resulting in a limit My > 611 GeV if B(b' — > Wt) = 1. For a b' quark 
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FIG. 5.8: The branching ratio of r' into the 
final states Wv, Zt, and hr, as a function 
of its mass, assuming Mh — 125 GeV. 




r 



200 300 400 500 600 



m t , [GeV] 



decaying only into Zb, there is an ATLAS search [68] based on 2.0 fb which results in My > 400 
GeV. In our case, we see from Figure [5761 that B(b' — > Zb) = 0.5 is a more likely scenario, in which 
case the limit from [(38| is about 360 GeV. However, the ATLAS analysis only uses Z — > e + e~ , so 
improvements can be expected both from using /i + /U~ and more integrated luminosity. As in the 
case of t'i, it would be useful to exploit the other decay modes in a comprehensive search strategy 
that allows the branching ratios to vary. In particular, the decay b' — > hb will lead to a nice signal in 
which there are at least 4 potentially taggable 6-jets. For example^ — > b'b' —> {Zb){hb) — > l + l~~bbbb 
should make for a background-free signal. 



The spectrum of the model we are considering also has an exotic lepton, the r', whose quantum 
numbers are those of a right-handed electron with its vector complement. If the r' decays promptly, 
it will be difficult to find. Assuming that mixing is only with the r, the branching ratios to final 
states Wv, Zt and hr are shown in Figure [5781 . The total width is determined by the e E coupling 
in eq. (|2.6p . but the branching ratios depend only on the r' mass. The production cross-section 
is rather low for this state, being electroweak strength, as is shown in Figure 15.91 However, one 
can produce unique signatures such as l + r~h + E™ 1SS that could be exploited at the LHC to 
simultaneously find the Higgs boson and the r' . A full exploration of these prospects will be 
pursued in another publication. 

If the t' is stable, it can be searched for as a weakly interacting heavy stable charged particle. 
The lifetime depends only on the mass and on the mixing coupling e E , with 



D. Search for r' at the LHC 




(5.11) 



We have taken the formal limit of M T i 3> Mh,Mz, M\y here for simplicity, and kinematic effects 
will lengthen cr by a factor of a few if the r' mass is not far above 100 GeV. Note that there is also 
an enhancement of the lifetime proportional to l/cos 2 /?, so that the r' could have a measurable 
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x' Mass (GeV) 

decay length with e E as large as a few times 10 -5 if tan j3 is large. While this may seem quite 
small, it is larger than the electron Yukawa coupling in the SM. 

The implications for the search are like that of a quasi-stable stau boson NLSP. The search 
strategies are very similar to that described in the t' section above, and so we shall not repeat it 
here. The result is that with 4.7 fb -1 of integrated luminosity, we can compare the cross-section vs. 
mass limits of [f33| to the cross-section computation in Figure 15. 9\ and conclude that there is a limit 
of quasi-stable r' mass of m T i > 450 GeV. We estimate that with 100 fb -1 of integrated luminosity 
at a 14TeV LHC phase, the reach for quasi-stable r' can extend up to nearly ITeV, which is well 
within the range of r' masses expected for Mh ~125 GeV assuming that Me = Mq = M\j at the 
unification scale, as illustrated by the examples of Figures 12.11 and 12.21 

VI. CONCLUSION 

A minimal GMSB model, with one 577(5) 5 + 5 messenger pair, can explain a Higgs mass of 
~125 GeV with even a sub-TeV gluino. This is accomplished by adding to the spectrum 10 + 10 
vector-like states, which then couple to the Higgs boson via the superpotential of eq. (|2.5p . The 
resulting radiative corrections can easily add 10 GeV or more to the light Higgs boson mass, which is 
crucial to achieve the ~125 GeV naturally, without requiring superpartners to be well above 1 TeV 
or invoking ad hoc non-GMSB stop mixing. We have paid special attention to cases inspired by 
unification of masses, Mq = Mjj, and mixing couplings, e v = e' v = e D , and we have characterized 
its parameter space. In this case, it is generic that the lightest exotic quark of the spectrum, t^, 
is mostly that with quantum numbers similar to right-handed top quark, with particular decay 
branching fractions. 

The most obvious implication for this scenario is the existence of low-scale super symmetry that 
should reveal itself at the LHC in the coming years. The searches for superpartners should follow 
the usual searches for GMSB models, which implies the existence of standard supersymmetry 
missing energy signatures with the addition of extra photons (or taus) if the NLSP is a neutralino 
(or stau) and decays promptly. The signals may feature also either the presence of the lightest 
Higgs boson h or a high multiplicity of taus due to wino decays in many events, depending on 
the messenger scale. If the decays of the NLSP are not prompt, the collider phenomenology will 
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be similar to that of standard scenarios with the neutralino being the LSP (i.e., stable NLSP on 
detector time scales), or there will be stable charged particle tracks from a quasi-stable charged 
NLSP stau. 

The scenario under consideration in this paper yields additional phenomenological implications 
due to the existence of the t\ 2 and b' and r' exotic fermion states. In previous sections we have 
explained that these states can also yield quasi-stable charged particle tracks, with sensitivity 
being nearly 1 TeV already for and nearly 0.5 TeV for r'. If the decays are prompt, the limits 
are reduced. In that case the t± pair-production signal is probably the most telling one for our 
scenario. We estimate sensitivity to the mass to be higher than 1800 GeV at 14 TeV LHC with 
more than 10fb _1 of integrated luminosity. If seen with the properties described in the previous 
sections, the signal would point to the existence of extra vector-like quarks that lift the Higgs boson 
mass to ~125 GeV. 



Appendix: Exotic quark and lepton couplings to W, Z, h and decay widths 

This Appendix is devoted to a systematic description of the interactions of quarks and leptons 
to the massive weak bosons W, Z, h, allowing for arbitrary flavor violation, and to formulas for the 
corresponding flavor- violating fermion decays. 

In the quark sector, we promote the third- family mixing parameters e U: e' v , and e D to couplings 
e^ 7 , ef, and ef respectively, where the index i = 1,2,3 indicates the three SM generations. The 
masses for up-type and down-type quarks in the gauge-eigenstate basis are then respectively 5x5 
and 4x4 matrices: 



VijVu qv u 
Mu k'v d 
\efv u kv u M Q J 



M d 



Vij v d 



Mr 



(A.l) 



where yfj and yfj are the 3x3 MSSM Yukawa couplings for the ordinary quarks, and the 
entries appear by a choice of basis. One can now obtain the gauge-eigenstate two-component left- 
handed fermions^ by applying unitary rotation matrices L, R, V and R' on the mass eigenstates 
Ui = (u, c, t, t±, t' 2 ) and Ui = (u,c,t,t l ,t 2 ), and eZj = (d,s, b, b') and Uj = (d,s,b,b ), so that 



L T M U R 
L' T M d R' 



diag(m M , m c , m t , m t / , m t > 2 ), 
di&g(m d ,m s ,m b ,m b ,). 



(A.2) 
(A.3) 



The first index of each of L, R, L' , R' is a gauge eigenstate index, and the second is a mass eigenstate 
index.-'' Then the interaction Lagrangian for couplings of W, Z, h to the quarks can be written as 



t We use the two-component fermion notations of . The four-component Dirac fields are 
* The notation used in [l2|] had a similar appearance but different index orderings. 
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(A.4) 



where the couplings for the boson are 
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(A.5) 
(A.6) 



and the couplings for the Z boson are 



d J d, 
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(A.7) 
(A.8) 
(A.9) 
(A.10) 



and the couplings for the lightest Higgs scalar boson are 



cos a(^L ki R nj y kn + L ki R^e v k + L 5i R nj e^' + L 5i R 4j k^ - sin a L Ai R 5j k' , (A.ll) 



1 



(A.12) 



The couplings to the heavier neutral Higgs bosons H° and A are obtained by the replacements 
(cos a, sin a) — > (sin a, — cos a) and (i cos /3, —i sin /?) respectively. 

Note that in the couplings of the W boson in eq. (|A.5|h the role of the SM CKM matrix is 
played by the restriction to the i,j = 1, 2, 3 subspace of the 5x4 matrix 



Kij — ~^^L* ki L' k j + L^L'^. (A. 13) 



k=l 



Clearly, neither the full matrix nor its restriction is unitary. (In the standard notation of 70], 
our K\\ is V u d, our K23 is V c b, etc.) Also, there is a nonzero coupling of the W boson to right- 
handed quarks in eq. (|A.6|h unlike in the SM. However, these flavor- violating effects do decouple 
as ef , ef and ef are taken to zero or as Mq and My are taken very large. Similarly, tree- level 
flavor-changing neutral currents of the Z boson couplings appear as the three terms with explicit 
reference to the exotic quarks' gauge-eigenstate indices 4,5 in eqs. (|A.7|) . (|A.8|) . and (jA.lOp . 

The widths of kinematically allowed flavor-changing two-body decays of quarks involving weak 
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bosons are given by 



T(u t ^Wd 3 ) = ^X 1 /\l,r w ,r j ){[l + r j -2r w + (l-r j )yr w ](\g^ d f + \g^ 



W |2\ 



32vr 
+ l2,/r. 



T(di -> Wuj) 



T(q t -> Z#) 



^A^(l )rwjri ){[l + r s - 2r w + (1 - r 3 f/r w ](\ g ^ + | 
+ l 2v ^Re[^4j}, 

^A^Cl.rz.rijjll + r, - 2r z + (1 - r^/rzldflft-.l 2 + lift, l') 

OZ7T L 9j 9j Hi 1j 

+12JT-Re[g Z i g%_ ]}, 

V J ^q\q^q\q^\' 

M, 



(A.14) 



(A.15) 



(A.16) 



with q = u or d, and A(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz, and rx = M x /Mg.. The special 
cases considered in the text above are t[ — > Wb, Zt, ht, and b' — > Wt, Zb, hb, both obtained by 
taking i = 4 and j = 3, with the mixing of exotic quarks to SM quarks restricted to the third 
family. The t[ decays were also discussed in [12j (using a different notation). 

In the case of a b' with My < My + Mw, there may be a competition between the two-body 
decays above and the flavor-preserving three-body decay through an off-shell W boson to SM 
fermions. In the approximation that flavor mixing between the exotic fermions and the SM leptons 
and first and second-family quarks is neglected, we obtain 



T(b' -> t'Jf 



M Q g 2 \V ff 
1536vr 3 



w „w 



where Vff is the standard CKM matrix for quarks (/ = u, c and /' = d, s) and is the Pontecorvo- 
Maki-Nakagawa-Sakata (PMNS) matrix for leptons (/ = neutrinos and /' = e, fi, r), and 



(i-^p 2 A 1 / 2 (l,g,^ i )A 1 / 2 (l,r / /x,r / Va) ^ 
iV rj +VTJ7) 2 [x - r w y + -yr w 



with 7 = r 2 ^ j Ml and r x = M\ /M§ , and 

fx = + r 4 / - x)X(x,r f ,r f >) + [(1 - r t / ) 2 - x 2 ][x 2 +x(r/ +r f ) - 2(r f - ?//) 2 ] j /a 



(A.19) 
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(A.20) 
(A.21) 



This formula is also valid (and smoothly approaches) the two-body decay width when the W boson 
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is on-shell, in the narrow-width approximation 7 -C rw 

1 



7T 



{x-r w ) 2 +^r w y/Ww 



5(x — ry/)- 



(A.22) 



The F\ kinematic part of this result was obtained in 71 ] . 
In the charged lepton sector, the 4x4 mass matrix is 



M E 



(A.23) 



where ef is a mixing coupling, with i,j = 1,2,3. The gauge eigenstate two-component fields are 
related by unitary rotations U,V acting on the mass eigenstate basis (e, /i,r, r') and (e, /J, r, t') in 
such a way that 



U 1 M. e V = diag(m e , m^, m T ,m 7 



(A.24) 



We assume that there are 3 light Majorana mass eigenstate neutrinos ^1,2,3, related to the gauge 
eigenstates v e » T by a unitary PMNS matrix N according to 
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The weak boson interactions with mass-eigenstate leptons are 



+ Z ^9 Z v kA^3 + 9j p .e\a% + gLM°%) + (y h e - ej h e t e J + c.c.), (A.26) 



where the couplings are 
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(A.28) 
(A.29) 
(A.30) 
(A.31) 



Note that unlike in the SM with 3 massive Majorana neutrinos, the effective PMNS matrix Mi 3 
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Y^k=i ^ki^kj is n °t unitary in general. The other change from the SM prediction comes from the 
left-handed coupling to the Z boson in eq. (|A.29h . This deviation from lepton universality is small 
in the limits that ef is small or Me is large. 

The resulting two-body decay widths for r' are 

r(r ' _> Wuj) = ^-(1 - r w )\2 + l/r w )\gWJ, (A.32) 
T(r' -+ Ze 3 ) = - r z f{2 + l/r*)|flf t | 2 , (A.33) 

OZ7T e 4 e J 

T(t> -> ^) = ^(1 - r fc ) a (|i^| a + |4, 4 | 2 ), (A.34) 

where the e 3 - = e,fi,T lepton mass is neglected for kinematic purposes, and the first decay should 
be summed over j = 1, 2, 3 when the neutrinos are not observed. In the numerical example in this 
paper and in [121], the special case is taken in which ej coupling is only non-zero for j = 3, so that 
electrons and muons do not mix with the r', and only the decays t' —> Wis, Zt, hr occur. 
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